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ABSTRACT 

Context. Following an extensive survey of the galactic plane by the INTEGRAL satellite, new hard X-ray sources are discovered with a signif- 
icant fraction of Cataclysmic Variables (CVs) among them. 

Aims. We report here the identification of one of these hard X-ray sources, IGR J00234+6141, as an accreting magnetic white dwarf of inter- 
mediate polar type. 

Methods. We analyse the high energy emission of the INTEGRAL source using all available data and provide complementary optical photo- 
metric and spectroscopic data obtained respectively in August and October 2006. 

Results. Based on a refined INTEGRAL position, we confirm the proposed optical identification. We clearly detect the presence of a 564 s 
periodic optical modulation that we identify as the rotation of the white dwarf. The analysis of the optical spectrum also demonstrates that the 
emission lines show a modulation in radial velocity with an orbital period of P^i-b = (4.033 ± 0.005) hr. 

Conclusions. The two periodicities indicate that IGR J00234-I-6I41 is a magnetic CV of the intermediate polar type. This is one of the faintest 
and hardest sources of this type detected by INTEGRAL. This confirms earlier conclusions that IPs contribute significantly to the population 
of galactic X-ray sources and represent a significant fraction of the high energy background. 
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1. Introduction 

The first complete galactic plane scan by the INTEGRAL 
satellite at energy above 20keV has shown that Cataclysmic 
Variables (CVs) constitute a significant fraction of the galactic 
high energy emission. The 2nd INTEGRAL source catalogue, 
from a ~10Ms dataset ( [Bird et aiT~||2006 ), lists 8 confirmed 
CVs out of a total of 209 detected sources, amounting there- 
fore to 4% of the catalogue. The last more sensitive survey, 
based on a 40 Ms exposure time (Bird et al. ||2007), has in- 



creased the catalog to 421 sources and more than doubled the 
total number of CVs (21), increasing their fraction to 5%. Most 
of the identified high energy CVs are Intermediate Polars (IPs) 



Send offprint requests to: J.M. Bonnet-Bidaud, email: bonnet- 
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* Based on observations obtained at the Haute-Provence 
Observatory (France) and at the Loiano Observatory (Italy) op- 
erated by the Istituto Nazionale di Astrofisica. 



i.e. accreting magnetic white dwarfs with intermediate (10^- 
lO^G) magnetic fields. Though the standard accretion column 
model predicts high energy emission of magnetic CVs (mCVs) 



(|Aizu [1973 I , only few sources have been studied above 20keV 



( Matt et al. [2000} [deMartino et al. |2001 1 [2004{ [Falanga et aL 
2005|. In the strongly magnetized (B> lO'G) polar systems, 
cyclotron cooling is an important mechanism to suppress the 
bremsstrahlung high temperature emission whilst it should be 
negligible among IPs. This could explain why most hard X-ray 
CVs are IPs. However there is so far no clear explanation of the 
high temperature achieved in some cases and this remains an 
open problem until a detailed characterization of each of these 
high energy CVs is obtained. Here we report on the identifica- 
tion of one of the very few CVs discovered by its high energy 
emission. 

IGR J00234-f6141 was first discovered from the INTEGRAL 



deep survey in the Cas-A region ( den Hartog et al. 2006 1. From 
a 1 .6 Ms exposure, a weak ( ~ 0.7 mCrab) source was detected 
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in the (20-50 keV) band with significant flux detected up to 
100 keV. The best source position was provided with an es- 
timated 3' accuracy, marginally consistent with a ROSAT soft 
X-ray counterpart, IRXS J002258.3-H6141 1 1, suggesting an X- 



ray binary ( den Hartog et al. 2006 1. Optical investigation of the 
much smaller («10") ROSAT error box led to the identifica- 
tion of a (R X! 16.5) object showing strong emission lines of 
the Balmer series typical of cataclysmi c variables (jH alpern & 
IMirabal |[2006l [Bikmaev et £171^^061 [Masetti et al. ||2006| ). 
From a preliminary short 2 hr photometric R-band observation, 
the presence of a possible 570 s periodic modulation was also 
reported suggesting a rotating white dwarf (Bikmaev et al. | 



2006 1. We present here the first detailed photometric and spec- 



troscopic observations leading to the secure identification of 
the source as an Intermediate polar. 

2. INTEGRAL analysis 



The INTEGRAL ( [Winkler et aL][2003l ) high energy emission 
(20-100 keV) of IGR J00234 -I-6141 was analyzed using all 
publicly-available IBIS/ISGRI ([ Ubertini et al.|2003|[Lebrun et 
|al.|2003| l data from March 2003 to October 2006. The data were 
collected for all ISGRI pointings within < 9° for a total effec- 
tive exposure ~ 3.6 Ms. To study the weak persistent X-ray 
emission, the ISGRI (20-100 keV) spectrum has been extracted 
from the mosaic images in four energy bands (20-28 keV), (28- 
40 keV), (40-57 keV) and (57-80 keV). The data reduction was 
performed using the standard Offline Science Analysis (OSA) 
software version 5.1. 

2.1. Source position 

The source is clearly detected in all energy bands with a max- 
imum significance of 6.8 cr in the (28-40 keV) band and a 
total (20-80 keV) significance of 8.9 cr. Based on the total 
exposure of ~ 3.6 Ms, the best position is determined to be 
a = 00''22'"54.3s, 6 = H-61°43'05" (J2000) with a 90% ac- 
curacy of 5.2'. This is consistent with the position derived in 
the third INTEGRAL catalogue from a lower 1.8 Ms exposure 
( [Bird et al. |2007) . The INTEGRAL position is now in better 
agreement with the proposed optical counterpart than the first 
determination by |den Hartog et al. ( 2006| l, reducing the separa- 
tion from 3.2' to 1.3' (seeFig.Q. 



We also note that contrary to what was reported by Bikmaev et 
aL 



(2006 



see their Fig. 6) the position of the ROSAT source, 
IRXS J002258.3H-614111, is fuUy consistent with the optical 
counterpart with an error circle much smaller than shown by the 
authors and an X-optical separation less than 6.1". The identi- 
fication of IGR J00234+6141with both the soft X-ray source 
and the optical star is therefore quite secure. 



2.2. Source spectrum 

Spectral analysis was performed with XSPEC version 12.0 
( Amaud|1996 Dorman et al. |200"3 i. A systematic error of 2% 
was applied to the ISGRI spectrum, which corresponds to the 
current uncertainties in the response matrix. All uncertainties 
in the spectral parameters are given at a 90% confidence level 




Fig. 1. New INTEGRAL position and uncertainty for the high 
energy source IGR J00234-I-6141, derived from our ~ 3.6 Ms 
exposure (cross and outer 5.2' circle). The inner (11" radius) 
circle is the ROSAT position, fully consistent with the proposed 
optical counterpart (arrow). 




Energy (keV) 



Fig. 2. The hard vFy energy spectrum of IGR J00234+6141 
with the best kT x 31 keV bremsstrahlung fit superposed and 
corresponding residuals (bottom). 



for single parameters. We first fit the spectrum using a sim- 
ple power-law model which gives a fit with;t'^/d.o.f. - 3.1/2, 
and the power-law index is determined to be F = 2.4^y^. A 
better fit was found with a bremsstrahlung model more com- 
monly used for CVs, this gave a ;^'^/d.o.f. - 2.6/2, with a 
bremsstrahlung temperature kT = 31.3^j4 keV. The flux in the 
(20-100 keV) range is 8.3 x 10"^^ erg cm^^ s"'. In Fig.|2] we 
show the unfolded spectrum and the residuals of the data to the 
bremsstrahlung model. 

We find significant diff'erences with the source parameters pub- 

( |2006l l. 



lished from a smaller dataset by [Barlow et al. | ( [2006| l, with 
a factor 2 higher bremsstrahlung temperature and with a factor 
1.3 lower flux. But our reduced ;t'^/d.o.f. = 2.6/2 is much more 
consistent than the one quoted by these authors (x^/d.o.f. - 
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6.1/3). For the energy band (20-40 keV), we also have a higher 
significance of 8.2cr (6.4cr for Barlow et al. 2006) that can be 
explained since we used only the pointings where the source 
was within < 9°, therefore increasing the S/N. 

3. Photometric observations 

IGR J00234-H6141 was observed on August 26 and 27, 2006 at 
the 1.5 m Loiano telescope (Bologna, Italy) equipped with the 



BFOSC detector operated with the Gunn g filter ( Gualandi & 

Merighi 1200111. The integration time was set to 30 s in all obser- 

vations. The CCD was windowed to a size of ~ 5 x 5 in order 

to reduce read-out time, giving a temporal resolution of 40 s. 

On August 26, the sky conditions were very poor, with winds 

ff 

and clouds and a seeing ~ 3 so that the source was observed 
for only ~40 min. On August 27, the sky conditions were pho- 
tometric with a seeing of 1.5-2 , although some cirrus were 
present during the night. We here then present the observations 
obtained on August 27. 

A total of 455 frames were acquired in 5.3 hr with the log of the 
observations reported in Table [T] The photometric data were 
reduced using /ra/ package with standard procedures includ- 
ing bias, flat-field and sky subtraction. Aperture 3" photometry 
was obtained for the target and several comparison stars. The 
flux ratios were obtained by dividing the counts of the target by 
the best combination (weighted mean) of three reference stars. 
These ratios were converted to fractional intensities by divid- 
ing by the mean flux ratio and heliocentric correction was also 
applied. The resulting light curve, shown in Fig.[3](top panel), 
displays clear short term periodic variations superposed on a 
longer term variability. A Discrete Fourier Transform (DFT) on 
the photometric time series shows a strong peak at 153 day"' 
(see Fig. [3] bottom panel). The low frequency part also shows 
some additional power at ~44 day"' and below 10 day"'. This 
latter could be due to a mixture of residual atmospheric extinc- 
tion and sky transparency fluctuations as well as possible con- 
tribution from variability at the orbital period 0.168 day (see 
below). We have simulated the long term variability using a 
third order polynomial function and two sinusoids at the fixed 
orbital frequency and another left free to vary. This latter is then 
found at43.88±0.23day-' (32.82+0.17 min). 
The DFT of the detrended light cui-ve is shown in the middle 
panel of Fig. |3] From a sinusoidal fit to the detrended curve, 
the main peak is found at 153.3 19±0. 167 day"'. A periodic 
563.53+0.62 s modulation is then clearly and safely detected 
with an ephemeris as: 



-pmax ^ HJD 2 453975.5 1782(7) + 0.006522(7) E, 



(1) 



where T™™ is the heliocentric anival time of the sine curve 
maximum. There is no indication of beat periods or harmonics. 
The light curve folded with this period is indeed very sinusoidal 
with an amplitude of (0. 100 + 0.003) mag (see Fig. |4]i. This pe- 
riodicity, typical of a spin period, indicates that the system is an 
intermediate polar. Though there appears to be an indication of 
a possible variability at 32.8 min, the identification of this type 
of longer period variations obviously requires a longer cover- 
age. 




0.5 0.55 
Time (HJD 2453975.0 + ) 



Spectral Window 




200 300 
Frequency (1/day) 

Fig. 3. Top : IGRJ00234+6141 g-band Ught curves (points) 
overlaid with a simulated curve (line) including the best spin 
period (564 s) and the long term variations (see text for de- 
tails). Middle ; Power spectrum (DFT) of the light curve, de- 
trended from the low-frequency variations. Ordinate are power 
amplitudes in units of 10"^. The inserted panel shows the spec- 
tral window Bottom : DTF of the observed light curve showing 
a strong peak at ~ 153 day"'. Some additional power is also 
present at lower frequencies (see text). 



4. Spectroscopic Observations 

IGR J00234+6141 was observed at the Observatoire de Haute 
Provence (France) in October 2006 (see Table [T]l. Long slit 
spectra were obtained with the Carelec spectrograph (Lemaitre 



|et al.|[T990| ), attached to the Cassegrain focus of the 193 cm 
telescope and using a CCD EEV (2048x1024 pixels) detector 
of 13.5 pixel size. A 133 A/mm grating was used, with a slit 
width of 2" , leading to a wavelength coverage of 3600-7200 A 
at a FWHM resolution of ~ 5.7 A. Exposure times were set at 
1128 s, twice the 564 s modulation, to smear any effect from 
the spin variability. The observations were obtained in nearly 
photometric conditions with only partial thin clouds during one 
night (21 Oct.). Seeing was typically of 2.5" to 3" . 

Standard reduction was performed in the ESO-MIDAS 
package, including cosmic rays removal, bias subtraction, flat- 
field coiTection and wavelength calibration. The wavelength 
calibration was checked on sky lines which were found within 
0.6 A from their expected wavelengths. All radial velocity 
measurements have been corrected from the Earth motion and 
from the small instrumental shifts measured on the OI 5577 A 
line and times have been converted in the heliocentric sys- 
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Table 1. Log of photometric and spectroscopic observations. 



Observations 


Range 


Res. 


Date 


UT(start) 


Exposure (min) 


Photometry 


g filter 


40 s 


2006 Aug. 27 


21:45 


320 


Spectroscopy 


3600-7200 


5.7 A 


2006 Oct. 21 


18:56 


244 




3600-7200 


5.7A 


2006 Oct. 24 


22:34 


134 




3600-7200 


5.7 A 


2006 Oct. 25 


18:24 


538 




Table 2. The parameters of the orbital modulation of the emis- 
sion lines radial velocities. 



563,5s Phase 



Fig. 4. The photometric g-band light curve folded with the best 
563.5 s ephemeris. The modulation is very sinusoidal with a 
full amplitude of 0. 10 mag. (full line). 




Frequency ( 1/DAY ) 

Fig. 5. Power spectrum of the Hq. radial velocities. The best 
value is found at Pn ~ 5.9 day"' (~4.1 hr). Secondary maxima 
are 1-day and 4-day harmonics. 



tem. Flux calibration has been performed using the standards 
BD+28 4211 and EG 247. 

The mean optical spectrum of IGR J00234+6141 is typi- 
cal of magnetic CVs with strong emission lines of the Balmer 
series, Hell (4686A), Hel (4471, 5875, 6678 and 7065A) and 
CIII-NIII (4655A), superimposed on a relatively blue con- 
tinuum. The mean V flux is estimated at ~ 7.5 x 10"'* 
ergs cm s A , which corresponds to a V magnitude of 



46.7, consistent to that reported by Bikmaev et al. (20061 



Lines 


y ( kms ') 


K(kms-') 


Phase (*) 




- 24.5 (5.2) 


47.1 (6.9) 


- 0.07 (0.07) 


HeII( 4686 A) 


-60.1 (8.7) 


71.3 (11.7) 


-1-0.08 (0.09) 




- 30.6 (4.5) 


53.9 (6.1) 


- 0.07 (0.06) 




-31.6 (3.0) 


60.0.(4.0) 


0.0 (fixed) 



* Phase of the blue-to-red zero crossing 
I-tr eiTor bars into parentheses 




Fig. 6. Radial velocities of the H^line folded with the orbital 
period Pn=4.033 hr. The best sine fit is also shown (Une). 



lines were measured using a single Gaussian least square fit 
procedure and were analysed to constrain the orbital period. 
Periodicities were searched using a Discrete Fourier Transform 
(DFT). The power spectrum for the strongest Hq. line gives the 
most accurate period determination (see Fig|5]l. The maximum 
power is seen at 5.94+0.14 day"' (4.03±0.10 hr) with sec- 
ondary peaks which are 1 d and 4d harmonics, corresponding 
to the night separations. The characteristics of the line modula- 
tion were further determined from a sine fit. The minimum 
corresponds to a period of Pn=4.033+0.005 hr, where the 
error bar is at a 99% level, computed for two independent pa- 
rameters. The orbital ephemeris for IGR J00234-I-6141 is de- 
termined as : 



Tn = HJD 2454033.4717(12) + 0.16804(21)E, 



(2) 



and [Masetti et al. | ( |2006[ l. The radial velocities for the main 



where Tn is the predicted heliocentric time of the blue-to-red 
radial velocity transition. 
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The radial velocities of the Hq, line, folded at the above 
best orbital ephemeris, are shown in Fig. [6] and the parameters 
of the corresponding best sine fits for the strongest emission 
lines are given in Table |2] It should be noted that a signifi- 
cant velocity scatter is present around ~ 0.25, at maximum 
red velocity. We carefully checked by using different sky lines 
that the effect is not due to a faulty wavelength calibration. In 
fact, the difference in velocity appears in the two consecutive 
cycles obtained in the same night (Oct. 25) when negligible 
wavelength shifts were observed. Inspection of the line profiles 
shows some distortions but an attempt to fit with two gaussians 
failed to give significant results. The B aimer lines show a ~ 
50-60 km s ' range in velocity amplitudes while the Hell lines 
display a slightly higher value with no significant phase shift. 

5. Discussion 

Out of the 21 sources identified as CVs in the 3rd INTEGRAL 
catalogue, IGR J00234+6141 is so far the unique confirmed IP 
that was discovered by its hard (E > 20 keV) X-ray emission. 
Apart for the close-by dwarf nova SS Cygni and the three po- 
lars (V834 Cen and the two asynchronous V1432Aql and BY 
Cam), all other sources are already known IPs with however 
five additional suspected CVs remaining to be identified from 
their period determinations (Bird et al. 2007 [ Barlow et al. 
2006} [Stephen et aL|[2006l ). Besides its INTEGRAL high en 
ergy emission, IGR J00234H-6141 has only one known low en 
ergy counterpart, the weak (0. 1-2.4 keV) ROSAT source, IRXS 
J002258.3+6141 1 1 ( jVoges et al.| 1999 ). From a search in the X- 
ray databases, a flux is also reported from the RXTE ALL Sky 
Monitor instrument, in the time interval from February 1996 to 
March 2007 ( |Levine et al.||1996] ). We extracted and analysed 
the (2-10 keV) light curve obtained at a 90 s resolution. The 
source is not significantly detected with a mean countrate of 
(0.050±1.65) c/s and a power spectrum of the light curve does 
not show any significant peak. 

Without any knowledge of the absorption, no attempt was 
made to fit joindy the ROSAT and INTEGRAL data. From 
the reported ROSAT hardness ratios, HR 1 = 1. 00+0.09 and 
HR2=0.51±0.16 ( |Voges etal.| 1999)1, a nd using the HR-Nh dis- 



tributions for CVs (see Motch et al. 1998| Fig. 4 and 5), a 



very rough estimate of the absorption can be obtained as Nh 
~ 2 X 10^' cm"-, significantly less than the total line-of-sight 
interstellar value (Nh ~ 7.3 x 10^' cm"^). With this Nh value, 
the extrapolation of the best INTEGRAL fit yields a (0.1-2.4 
keV) flux of 8.3 X 10"'^ ergs cm"^s"^(with nearly a factor 2 
uncertainty) . Using the count conversion given in Voges et al. 
1999) 1, the estimated (0.1-2.4 keV) ROSAT flux is significantly 



lower at (1.1 ± 0.2) x 10~ ergs cm"-s" which may indicate 
that a more significant absorption is present. 
IGRJ00234+6141 is up to now the faintest CV detected 
by INTEGRAL. There is no good indication of the distance 
though a useful lower limit can be set using the infrared 
magnitudes from the 2MASS survey and the CV donors se- 
quence recently computed by Knigge ( 2006[ ). For an orbital 
period of 4.03 hr, the donor is expected to be an M3 star 
with an absolute magnitude of Mk ~6.2. If attributed solely 
to the secondary, the observed K magnitude (K=14.77+0.12) 



( jSkrutskie et al.|)2006) l would then yield a minimum distance 
of ~520 pc. To compare with the other high energy IPs de- 
tected by INTEGRAL, a distance independent colour-colour 
diagram was built by comparing the flux ratios in the ROSAT 
(L=0.1-2.4keV), and the low (M=20-40keV) and high (H=40- 
100 keV) INTEGRAL energy range s (see Fig. [7 i. Flux ra- 



tios are derived from the ROSAT (Voges et al. 



T999) l and 



INTEGRAL ( jBird et al. |2007) l catalogues for uniformity rea- 
sons though fluxes derived from the actual spectral fits may 
sHghtly diff'er for IGR J00234+6 141 (this analysis), V709 Cas 
(Falanga et al. 2005) and IGR J15479-4529 (|Bonnet-Bidaud 



[& Falanga ||2007( ) for which a more detailed analysis exists. 
The sources shown in Fig. [7] are the confirmed IPs with known 
spin and orbital periods. With H-M= (0.58 ± 0.17) and M-L= 
(7.23 + 2.40), IGR J00234+6141 has characteristics compara- 
ble to the hardest sources. 

The high observed temperature (kT « 31 keV) is a lower limit 
to the maximum shock temperature and therefore indicative 
of a white dwaif (WD) mass higher than ~ 0.75Mq (see 
|Aizu |1973) l. Assuming that the Balmer lines trace the WD or- 
bital motion, the radial velocity amplitude gives a system mass 
function of f„, = (3.7+0.7) xlO"^ M©. Together with the WD 
minimum mass, this yields an estimate of the system inclina- 
tion i > (21 - 27)° for a companion mass in the range (0.3- 
0.7)Mq comptatible with the orbit and a maximum inclina- 
tion of 37° for a 1.4 M© WD mass. If a distance of ~ 500 
pc is assumed, the unabsorbed luminosity will be L ,(0.1 -100 
keV)~ 6.7 X 10^^ (d/500pc)^ ergs s"', corresponding to a rather 
low accretion rate of M ~ 0.51xlO"'g s^'fora- O.75M0 WD. 
IGR J00234+6141 , with a ratio P„/Pa=0.04, is also in a rela- 
tively high degree of asynchronism. For the above WD mass 
and accretion rate, the condition for accretion to take place, 
corresponding to a magnetospheric radius lower than the coro- 
tation one, will imply a magnetic moment lower than 1 .0 x 10^^ 
Gcm^, corresponding to a surface magnetic field of ~ 2.4 x 10^ 
G. Unless the white dwarf is significantly less massive than 
indicated by the shock temperature, IGR J00234+6141 could 
therefore be a particular case of a low accreting IP with a low 
magnetic field. 



6. Conclusions 

The analysis of INTEGRAL data has shown that 
IGR00234+6141 is among the hardest X-ray emitting 
CVs with a probable low magnetic field. We have identified 
an orbital variability in the radial velocities of emission lines 
at a 4.033 hr period. We have also found a pulsation with a 
period of 563.5 s and an optical amplitude of 0.10 mag which 
is a clear sign of accretion onto a magnetic white dwarf. This 
period is identified as the rotational period of the compact 
object and should reveal itself as strong X-ray pulsations. The 
pure sinusoidal modulation and the absence of any harmonics 
is an indication of a one-pole accretion. Due to the faintness of 
the source, the INTEGRAL data do not allow a detailed study 
of the variability so that further X-ray observations are clearly 
needed to confirm these characteristics. 
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V1223Sgr RX2133 



M-L [20-40]/[0. 1-2,4] keV 

Fig. 7. High-energy colour-colour diagram of the confirmed IPs 
detected by INTEGRAL. The energy bands are indicated in the 
labels. Sources for which only upper values exist are marked by 
crosses. IGR J00234-H6141(filled square) is among the hardest 
sources. Note that FO Aqr with M-L= 59.8, due to significant 
absorption, is kept out of the scale. 
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